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ABSTRACT 

We present the first comprehensive analysis of Hubble Space Telescope (HST) observations of short-duration 
gamma-ray burst (GRB) host galaxies. These observations allow us to characterize the galactic and local 
environments of short GRBs as a powerful constraint on the nature of their progenitors. Using the HST data 
for 10 short GRB hosts we determine the host morphological properties, measure precise physical and host- 
normalized offsets relative to the galaxy centers, and study the locations of short GRBs relative to their host 
light distributions. We find that most short GRB hosts have exponential disk profiles, characteristic of late-type 
galaxies, but with a median size that is twice as large as that of long GRB hosts, commensurate with their 
higher luminosities. The observed distribution of projected physical offsets, supplemented by ground-based 
measurements, has a median of rj 5 kpc, about 5 times larger than for long GRBs, and in good agreement with 
predicted offset distributions for NS-NS binary mergers. For the short GRB population as a whole we find 
the following robust constraints: (i) > 25% have projected offsets of < 10 kpc; and (ii) > 5% have projected 
offsets of > 20 kpc. We find no clear systematic trends for the offset distribution of short GRBs with and 
without extended soft emission. While the physical offsets are larger than for long GRBs, the distribution of 
host-normalized offsets is nearly identical due to the larger size of short GRB hosts. Finally, unlike long GRBs, 
which are concentrated in the brightest regions of their host galaxies, short GRBs appear to under-represent the 
light distribution of their hosts; this is true even in comparison to core-collapse and Type la supernovae. Based 
on these results, we conclude that short GRBs are consistent with a progenitor population of NS-NS binaries, 
but partial contribution from prompt or delayed magnetar formation is also consistent with the data. Our study 
underscores the importance of future HST observations of the larger existing and growing sample of short GRB 
hosts, which will allow us to delineate the properties of the progenitor population. 
Subject headings: gamma-rays:bursts 



1. INTRODUCTION 

The galactic and local environments of cosmic explosions 
provide powerful insight into the nature of their progenitors. 
For example, past studies of supernova (SN) environments 
have demonstrated that Type la and Type Ib/Ic/II events arise 
from distinct progenitor systems since the former are located 
in all types of galaxies, while the latter occur only in star 
formi ng galaxies, pointing to a d irect link with massive stars 
(e.g., Ivan den Bergh et all 120051) . In a similar vein, long- 
duration gamma-ray bursts (GRBs; duration, T90 > 2 s) have 
been linked with massive stars through t heir exclusive asso- 
ciation with star form i ng galaxies (e g.. IBloom et al] 119981 : 
iDjorgovski et alj|1998t iFruchter et all 1 1 9991) . Short-duration 
GRBs, on t he other hand, are now known to res i de in all types 
of galaxies feerger et alJ | 2QQ^|FjaxsLal] |2005l: iGehrels etall 
12005b IBloom et al.ll2006t lBergerll2009h . Moreover, even the 
star forming host galaxies of short GRBs differ from those of 
long GRBs; they have higher luminositi es and metall icities, 
and lower specific star formation rates (Berger 2009). The 
difference between long and short GRB host ga laxies, along 
with the lack of SN associations for short G RBs (Berg eret alJ 
200HlFox et alJl200llGehrels et alJl200"llHjorth et alJl2005b 
Bloom etaU l2006t ISoderberg et al.l 120061) . demonstrate that 
they have distinct progenitor populations. In particular, at 
least some short GRBs are associated with an older progenitor 
population. 
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Equally important are the local, sub-galactic environments. 
In the case of long GRBs, the distribution of their projected 
physical and host-normalized offsets relative to the host cen- 
ters matches the overall exp ected distribution for massive 
stars in an exponential disk (Blo om et al.l 12002). An anal- 
ysis of the brightness distribution at the locations of long 
GRBs indicates that they are disproportionately concentrated 
on the brightest regions of their hosts, primarily in compari- 
son to core-collapse SNe , which follow the o verall light dis- 
tribution of their hosts (IFruchter et al.l l2006). Both of these 
studies have relied on high angular resolution Hubble Space 
Telescope (HST) imagin g observations. PreUminary stud- 
ies of short GRB offsets dBerger et aT]l2005t iFox et alJl2005b 
Bloom & Proch aska 2006; Soderberg et al. 2006; iTroja et all 
2008tlD , Avanzo et al.l2009l) reveal somewhat larger projected 
physical offsets than for long GRBs, and have also led to a 
claimed trend of smaller offsets fo r short GRBs with extended 
X-ray emission (Tro ia et al]|2008l) . No study of the locations 
of short GRBs relative to their hosts light distribution has been 
published so far. 

Progenitor models of short GRBs lead to distinct predic- 
tions about the distribution of host properties and the local en- 
vironments of short GRBs (measured by their offsets and lo- 
cation relative to the host light distribution). In particular, the 
popular model of ne utron star and/or black hole binary merg- 
ers (NS-NS/NS-BH: lEichler et Ml9&k iNaravan et al.lfT992l) 
predicts larger offsets than for the massive star progenitors of 
long GRBs due to potential systemic velocity kicks. Various 
authors have employed population synthesis models to predict 
the distribution of offsets by convolving distribut ions of kick 
velocities, merger timescale, and galaxy masses (IBloom et alJ 
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IT991 iFrver et1d1[l99l iBelczvnski et al.ll2006h . Fo r Milky 
Way mass galaxies, appropriate for short GRB hosts (Berger 
2009), the predicted offset distributions have a median of 
~ 5 — 10 kpc, with a broad tail extending to tens of kpc. On 
the other hand, progenitor models that invoke magnetars, ei- 
ther from a young population or through delayed formation 
in a WD- WD merger or white dwarf accretio n-induced col- 
lapse (iLevan et al.l l2~006b; Me tzger et~ai1l2008l) . are expected 
to have a modest offset distribution since these systems do not 
experience kicks. Similarly, any systematic differences in the 
progenitors of short GRBs with and without extended X-ray 
emission may be revealed in the offset distribution and spe- 
cific sub-galactic environments. 

To confront these models with observations we require high 
angular resolution imaging, best provided by HST. Such ob- 
servations provide detailed information on the host galaxy 
morphological properties (e.g., exponential disk versus de Va- 
coulers profile, effective radius), as well as the ability to pre- 
cisely measure offsets and the distribution of short GRBs rel- 
ative to their host light. HST observations have served as the 
backbone for detailed studi es of long GRB environments and 
host galaxy morphologies dBloom et aT]|2002l; iFruchter et al.1 



1200a Wain wright et al. 2007). To date, HST obser vations of 
only 3 short GR Bs have been publis hed (05070 9: iFox et al 



2005; 060121: ILevan et all 120061 080503: iPerlev etal 



2009), and in only one case (050709) was the host morphol- 
ogy addressed. 

Here we present the first comprehensive analysis of all 
short GRB host galaxies observed with HST to date 3 . Us- 
ing these observations we determine the host morphologies 
and structural properties ($3), we calculate precise physical 
and host-normalized offsets using accurate astrometry relative 
to ground-based afterglow images (®, and we construct the 
first distribution of short GRB locations relative to their host 
light (©. We draw conclusions in the context of progenitor 
models in $6] Throughout the paper we compare and contrast 
the results of our analysis with similar studies of long GRBs. 

We find that: (i) the short GRB hosts have systematically 
larger effective radii than long GRB hosts, in good agreement 
with their higher luminosities; (ii) the observed short GRB 
projected physical offset distribution has a median of about 5 
kpc, about a factor of 5 times larger than long GRBs, while 
for the population as a whole > 25% have projected offsets of 
< 10 kpc, and > 5% have projected offsets of > 20 kpc; (iii) 
both the observed physical offset distribution and the robust 
constraints closely match the predicted offset distribution of 
NS-NS binaries; and (iv) short GRBs uniformly trace the light 
distribution of their hosts, similar to core-collapse SNe, but 
distinct from long GRBs. 

2. HUBBLE SPACE TELESCOPE DATA REDUCTION AND ANALYSIS 
2.1. Sample 

We present HST optical observations of ten short GRB host 
galaxies obtained with the Advanced Camera for Surveys 
(ACS) and the Wide-Field Planetary Camera 2 (WFPC2). The 
data were obtained as part of programs 10119, 10624, and 
10917 (PI: Fox), as well as 10780 and 11176 (PI: Fruchter). 
These programs targeted all short GRBs with optical and 
X-ray (XRT) positions from May 2005 to December 2006, 

3 We do not repeat the analysis for GRB 080503 since no convincing host 
galaxy was identified jPerlev et al.ll2009l) . Future HST ACS or WFPC3 ob- 
servations will provide better constraints on an underlying host than available 
from the existing WFPC2 observations. 



which were visible during HST 2-gyro operations. In this 
time frame, the only short burst that was not observed was 
GRB 060801. Thus, the sample in this paper is nearly com- 
plete in relation to the short GRBs with optical/X-ray af- 
terglows 4 . The HST observations of G RBs 050709 and 
060 121 have been pu blished previously by IFox et alJ (120051) 
and lLevan et al.l (l2006al) . respectively, but we re-analyze them 
here in a uniform fashion and perform a more comprehen- 
sive analysis of the host morphology and burst environment. 
Seven of the nine bursts have been localized to sub-arcsecond 
precision from optical afterglow detections, and five have 
known redshi fts (Table Q}; we us e a constraint of z > 1 .4 for 
GRB 051210 (IBerger etalJliOOlh . 

Six of the seven short GRBs with sub-arcsecond localiza- 
tion are robustly associated with host galaxies 5 (we present 
a host identification for GRB 060313 in this paper; see Ap- 
pendix [A}. The sole exception is GRB 061201, for which 
we explore two possible host galaxy associations based on 
the HST observations (see Appendix EJ; previously only 
one galaxy ( at z = 0.111) was considered a potential host 
( IBerger et all 120071; IStratta et al.ll2007l) . Details of the GRB 
properties and the HST observations are provided in Table Q] 

Throughout the paper we use the standard cosmological pa- 
rameters, Hq = 71 km s" 1 Mpc" 1 , f2 m = 0.27, and S7a = 0.73. 
All r eported magnitudes ar e corrected for Galactic extinction 
using Schlegel et al. ( 1998) dust maps. 

2.2. Hubble Space Telescope Data Reduction 

We retrieved pre-processed images from the HST archive 6 
for all available short GRBs. We distortion-corrected and 
combined the ind ividual exposures using the IRAF task 
mult idrizzle (IFruchter & Hooldl2002n koekemo er et al.1 
2002). For the ACS images we used pixfrac = 0.8 and pixs- 
cale = 0.05 arcsec pixel" 1 , while for the WFPC2 images, we 
used pixfrac=1.0 and pixscale = 0.0498 arcsec pixel" 1 , half 
of the native pixel scale. The final drizzled images, flux- 
calibrated to the AB magnitude system according to the ACS 
and WFPC2 zeropoints, are shown in Figures [TlfTOl 

2.3. Surface Brightness Profile Fitting 

We use two methods to fit the surface brightness profiles 
of the short GRB host galaxies. F irst, we use the gal fit 
software package (Peng et al. 2007) to construct the best two- 
dimensional ellipsoidal model of each galaxy image. Second, 
we use the IRAF task ellipse to produce elliptical inten- 
sity isophotes and to construct one-dimensional radial surface 
brightness profiles. We further use ellipse to measure the 
integrated AB magnitude for each galaxy (listed in Table[T|). 

2.3.1. galfit 

As an input to galfit we generate point-spread-function 
(PSF) models for each instrument and filter combination us- 
ing the tinytim software package. We assume a constant 
spectrum in F„; the difference in the 90% encircled energy 

4 We do not include in this analysis GRBs 051227, 060505, and 060614 
all of which had durations well beyond 2 s, even when isolating the initial 
"short" emis sion episode. The HS T data for GRBs 060505 and 060614 were 
published in lOfek et ail 120071) and lGal-Yam et all 420071) . respectively. 

5 For a complete discussion of the galaxy associat i ons we refer the reader 
to the following p apers: GRB 050 709: [Fox et all (200% GRB 050724: 
I Berger et alJ f2005): GRB 051221: Soderberg et al] (2009); GRB 060121: 
lLevanetalJ <2006a): GRB 061006: IBerger et all <2007[) and fD'Avanzo et all 
120091) . 

6 see http://archive.stsci.edu/hst/ 
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width of the PSF for a spectral index ranging from -2 to is 
only 1%. We additionally correct for distortion in the ACS 
instrument, and use a sub-sampling factor of 2 for WFPC2 to 
appropriately account for the reduced pixel scale in the driz- 
zled images. 

In each observation we fit the host galaxy image with a PSF- 
convolved Sersic profile: 

S(r)=E c exp{- K „[(r/r f ,) 1 /"-l]}, (1) 

where n is the concentration parameter (n = 1 is equivalent to 
an exponential disk profile, while n = 4 is the de Vaucouleurs 
profile), k„ ks In- 1/3 + 4 /405 w + 46/255 15rc 2 is a co nstant 
that is coupled to the value of n dCiotti & Bertin| [l999)). r e is 
the effective radius, and E c is the effective surface brightness 
in flux units. In the subsequent discussion, tables, and figures 
we use surface brightness in units of mag arcsec" 2 , designated 
as \i e . 

In all cases we fit the host galaxies with a single 7 Sersic 
profile and allow the parameters to vary freely. The resulting 
best-fit values of n, r e , and /i e , as well as the integrated AB 
magnitudes are provided in Table [2] For host galaxies that 
are detected at a low signal-to-noise ratio we find that a wide 
range of n values can account for the observed morphology. 
In these cases we fit the host galaxies with n fixed at values of 
1 and 4, and provide the results of both models in Table [2] 

The galfit models and residual images for all instru- 
ment/filter combinations are shown in Figures QJ9] Objects 
for which both n = 1 and n = 4 models provide an adequate fit 
are shown for both cases. 

2.3.2. Radial Profiles from /^MF/ellipse 

We use ellipse to generate elliptical isophotes for each 
host galaxy, with the center and ellipticity of each isophote 
allowed to vary 8 . The resulting radial surface brightness pro- 
files in units of AB mag arcsec" 2 are shown in FigureQT] We 
fit each profile with a Sersic model (Equation [TJ using n, r e , 
and jit e as free parameters. The best-fit values are listed in 
Table [2] and the resulting models are shown in Figure QT| 
We find adequate fits in all cases, although some host galax- 
ies clearly exhibit radial complexity due to irregular structure 
and/or an edge-on orientation. 

2.4. Astrometry 

To determine the location of each short GRB relative to 
its host galaxy we perform differential astrometry using op- 
tical and near-IR images of the afterglows 9 . With the excep- 
tion of GRB 050709, whose afterglow is directly detected in 
HST/ACS observations, we use ground-based images from 
Magellan, Gemini, and the VLT. The astrometric tie be- 
tween the afterglow and host images is performed using point 
sources in common between the two images; the source of the 
afterglow image and the number of astrometric tie objects are 
listed in Table[3] In the case of ground-based to HST astrom- 
etry we use a range of 15 - 85 common objects, with the num- 
ber depending on the density of stellar sources in the field, the 
depth of the images, and the field-of-view. To determine the 

7 In the case of GRBs 050509b, 050709. 060121, and 061006 we 
fit additional Sersic and point-source components to account for fore- 
ground/background objects. These components are not considered to be part 
of the host galaxy. 

8 For the hosts of GRBs 060121 and 060313, with low signal-to-noise 
detections, the center and ellipticity were held fixed throughout the fit. 

' Optical afterglows have not been detected in the case of GRBs 050509b 
and 051210. 



astrometric tie we use the IRAF astrometry routine ccmap. 
We find that a second-order polynomial, with six free param- 
eters corresponding to a shift, scale, and rotation in each co- 
ordinate, provides a robust astrometric tie in all cases. The 
resulting rms values are ctgb^hst = 13-30 mas (Tabled. 

We next consider the uncertainty in the afterglow position 
from each ground-based image. The centroiding accuracy 
depends on the size of the PSF and the signal-to-noise ratio 
(S/N) of the afterglow detection, ag = #fwhm/(S/N). We de- 
termine <76),grb for each GRB using the SExtractor pro- 
gram 10 (Table 0). In the case of GRBs 050724 and 051221 
we find that <70,grb is significantly smaller than <7gb->hst; 
for GRBs 060121, 060313, and 061006 the two sources 
of uncertainty are comparable; and for GRB 061201 ctsgrb 
dominates. The afterglows of GRBs 050509b, 051210, and 
060502b have only been detected in X-rays, with the Swift 
X-ray Telescope (XRT), and as a result their positional uncer- 
tainty is O-0.GRB ~ 1.7-5.8" ( Table E} . We n ote t hat the XRT 
positio ns from the catalogs of Butlei] (120071) and lEvans et all 
(2009) exhibit relative offsets of 5.6" (GRB 050509b), 3.0" 
(GRB 051210), and 2.9" (GRB 060502b) suggesting that the 
true positional uncertainties (including systematics) are larger 
than their quoted statistical uncertainties. 

The final source of uncertainty in the relative position of the 
GRB and host galaxy is the centroiding accuracy of the host in 
the HST images. To determine this uncertainty we again use 
SExtractor. The resulting values of ag ga i are listed in Ta- 
ble[3] We find that for GRBs 050724, 051221, and 061006 the 
host centroid uncertainty is smaller than both ogb^hst and 
ce.GRB, while for GRBs 060121 and 060313 cre,gai is compa- 
rable to ag grb- The combined offset uncertainties are listed 
in Table [3] ' 

A much more accurate relative position is available for 
GRB 0507 09 since the a fterglow was detected in HST/ACS 
images dFox et al.l 120 05). The limiting factor is thus ag for 
the afterglow and host. We find from the first HST/ACS ob- 
servation (2005 July 15.6 UT) that ag. grb = 1 mas, while 
C0,gai = 1.4 mas. In addition, we also astrometrically tie the 
final epoch of ACS imaging and the WFPC2 image to the first 
epoch. Since the afterglow is no longer detected in these im- 
ages, this allows us to study the burst environment. For the 
final ACS image we find cthst^hst = 8 mas, while for the 
WFPC2 image we find cthst^hst = 14.5 mas. These uncer- 
tainties clearly dominate over the centroiding errors of the 
afterglow and host galaxy. We do note, however, that the 
complex morphology of the host galaxy (@ introduces a sys- 
tematic uncertainty in the definition of the host "center". By 
varying the signal-to-noise threshold in SExtractor, we 
find that the centroid of the host shifts by as much as 50 mas, 
which dominates over the statistical uncertainty in the source 
position. 

Similarly, we find a more accurate offset for GRB 050724 
from a detectio n of the afterglow and host in ground-based 
near-IR images (Ber ger et al.l l2005). The combined afterglow 
and host centroid uncertainty in these images is about 6 mas, 
compared to a total uncertainty of 20 - 25 mas for the HST 
images. 

2.5. Host Light Distribution 

To determine the brightness of the GRB location relative 
to the host lig h t distr ibutio n, we follow the methodology of 
iFruchter et all (120061) and iKellv et al.l (120081) and calculate 

10 |http: //sextractor. sourcef orge ■ net/| 
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from each galaxy image the fraction of total light in pixels 
fainter than the afterglow position. Six bursts have differential 
astrometric positions of better than 1 pixel (050709, 050724, 
051221a, 060121, 060313 and 061006; Table EJ. For each 
image we create an intensity histogram of a 5" x 5" region 
centered on the host galaxy and determine a lu cut-off level 
for the host by fitting a Gaussian profile to the sky brightness 
distribution (equivalent to a signal-to-noise ratio cut-off of 1). 
We then plot the pixel flux distribution above the appropriate 
cut-off level for a region surrounding the host, and determine 
the fraction of light in pixels fainter than the afterglow pixel; 
see Table |U 

3. MORPHOLOGICAL ANALYSIS 

Using the results of the gal f it analysis and the radial sur- 
face brightness profiles we first classify the short GRB hosts 
in terms of their Sersic n values. From the gal fit analysis 
we find that three hosts (GRBs 050709, 051221a, and 061006) 
are best modeled with n ~ 1, corresponding to an exponential 
disk profile, while two hosts (GRBs 050509b and 050724) 
are best modeled with n w 3 and w 5.6, respectively, typi- 
cal of elliptical galaxies. We note that GRB 050724 possibly 
exhibits weak spiral structure, whi ch may explain the re sult- 
ing value of n w 3 (see Figure[3]and Malesani et alj200~7l) . but 
this putative spiral structure is clearly sub-dominant relative to 
the elliptical structure. The final three hosts (GRBs 051210, 
060121, and 060313) are equally well modeled in gal fit 
with a wide range of n values, and we provide results for both 
n = 1 and n = 4 in Table[2] 

We find identical results using Sersic model fits to the one- 
dimensional radial surface brightness profiles generated with 
ellipse (Figure [TT| and Table |2j. However, with this ap- 
proach we find best-fit values of n ~ 1 for the three host galax- 
ies with ambiguous gal fit results, suggesting that they are 
indeed better modeled as exponential disks. We therefore con- 
clude that of the eight short GRB host galaxies studied here 
only two can be robustly classified as elliptical galaxies based 
on their morphology. A similar fraction was determined in - 
dependently from spectroscopic observations (Berger 2009). 
The distribution of n values is shown in Figure fT2l 

As can be seen from the gal fit results, the Sersic models 
of the two elliptical hosts exhibit significant residuals (Fig- 
ure sQ] and [3]). This is a well-known effect for bright elliptical 
galaxies, which generally require a multi-parameter power- 
law plus Sersic fit that account s for a flatter core th an expected 
in the de Vaucouleurs model dTruiillo et al. 2004). Since we 
are here mainly interested in the distribution of n values and 
a comparison to long GRB hosts, we retain the simple Sersic 
formulation. 

We also find significant residuals for a one-component Ser- 
sic fit of the host galaxy of GRB 050709, which has an irreg- 
ular morphology dominated by an exponential profile (Fig- 
ure |2). This is the only clearly irregular galaxy in the sam- 
ple. Finally, we find that the hosts of GRBs 051210, 060121, 
and 061006 exhibit significant bulges, clearly seen in their ra- 
dial surface brightness profiles (Figure ITTb. For the host of 
GRB 061006, which was observed in two filters, the bulge 
component is more significant in the F814W filter than in 
the F555W filter, as expected for an older stellar population; 
the burst appears to coincide with this bulge component (Fig- 
ured. 

The gal fit and radial profile fits also yield values of the 
effective radius, r e , for each host galaxy. We find a range 



of w 0.2-5.8", corresponding to physical scales 11 of about 
1.4-21 kpc. The smallest effective radius is measured for 
the host of GRB 060313, while the host of GRB 050509b has 
the largest effective radius. The median value is r e m 3.5 kpc. 
We adopt the best-fit values from the radial surface brightness 
profiles, and plot the resulting distribution, as well as r e as a 
function of n, in Figure [T2l 

Finally, the effective surface brightness values range from 
/j, e w 21 to « 27 AB mag arcsec" 2 . The galaxy with 
the highest surface brightness is the host of GRB 050724, 
while the lowest surface brightness is measured for the 
host of GRB 060121. The integrated magnitudes range 
from about 16.3 AB mag (GRB 050509b) to 26.4 AB mag 
(GRB 060313). 

3.1. Comparison to Long GRB Host Galaxies 

A comprehensive morphological analysis of long GRB host 
galaxies using HST observations with the STIS, WFPC2, 
and AC S in struments has been carrie d out by lConselice et al. 
(2005) and IWainwright et al.l d2007l) . In Figure [12] we com- 
pare the values of n and r e measured for long GRB hosts by 
IWainwright et alJ (120071) to the values measured here for short 
GRB hosts. Two clear trends emerge from this comparison. 

First, all long GRB hosts have n < 2.5, and the median value 

i n 1 — — 1 —i 

for the population is (n)»l.l (IWainwright et al. 2007). Thus, 
all long GRB hosts are morphologically classified as expo- 
nential disks, while 2 of the 8 short GRB hosts studied here 
exhibit de Vaucouleurs profiles. However, for the hosts with 
n < 2, the distributions of n values for both populations appear 
to be similar. 

Second, short GRB hosts have larger effective radii, with 
(r e ) k 3.5 kpc, compared to (r e ) w 1.7 kpc for long GRB 
hosts ( Wainwri ght et alJ l2007). A Kolmogorov-Smirnov (K- 
S) test indicates that the probability that the short and long 
GRB hosts are drawn from the same underlying distribution 
of host galaxy effective radii is only 0.04. If we remove from 
the sample GRBs 050509b and 05 1 2 1 (which have only XRT 
positions) we find that the K-S probability is still only 0.09. 
Thus, we conclude with high significance that short GRB host 
galaxies are systematically larger than long GRB hosts, and 
that this result is not affected by host associations based on 
XRT positions. The larger sizes of short GRB hosts are ex- 
pecte d in the contex t of the galaxy size-luminosity relation 
(e.g.. lFreemarifl970l) . We recently showed that short GRB 
hosts are systematically mo re luminous t han long GRB hosts 
by about AMb ~ 1.2 mag (Berger 2009) and therefore their 
sizes are expected to be correspondingly larger. 

An additional striking difference between the hosts of long 
and short GRBs is the apparent dearth of interacting or ir- 
regular galaxies in the short GRB sample. Of the eight 
host galaxies studied here, we find only one irregular galaxy 
(GRB 050709) and none that appear to be undergoing merg- 
ers. In contrast, the fraction of long GRB hosts with an ir- 
regular or merger/interac tion morphology is about 30-60% 
(IWainwright et alJ l2007h . The interpretation for this high 
merger/interaction fraction in the long GRB sample is that 
such galaxies are likely undergoing intense star formation 
activity triggered by the merger/interaction process, and are 
therefore suitable sites for the production of massive stars. 

11 For the faint hosts without a known reds hift (GRBs 051210 , 060121, 
060313, and possibly 061201) we assume z = 1 I Berger et al. 2007), and take 
advantage of the relative flatness of the angular diameter distance as a func- 
tion of redshift beyond z ~ 0.5. 
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The lack of morphological merger signatures in the short 
GRB sample indicates that if any of the hosts have undergone 
significant mergers in the past, the delay time between the 
merger and the productio n of a short GRB is > 10 9 yr (e.g., 
iBarnes & Hernquistlll992l) . 

4. OFFSETS 

We next turn to an analysis of short GRB offsets relative to 
the centers of their host galaxies. Based on the astrometric tie 
of the HST host observations to ground-based afterglow ob- 
servations, we find that the projected offsets are in the range of 
» 0.12- 17.7" (Tabled. The corresponding projected physi- 
cal offsets are about 1-64 kpc, with a median value of about 
3 kpc. The largest offsets are measured for GRBs 050509b 
and 051210, but these are based on Swift/XKT positions with 
statistical uncertainties of about 12 and 18 kpc, respectively 
( and possibly larger if we consider systematic uncertainties; 
{32.41 i. If we consider only the bursts with sub-arcsecond af- 
terglow positions we find that the largest offset is 3.7 kpc 
(GRB 050709), and that the median offset for the 6 bursts 
is 2.2 kpc. In the case of GRB 061201 the host association 
remains ambiguous (see Appendix IB1. but even for the near- 
est detected galaxy the offset is about 14.2 kpc. The obvious 
caveat is that an undetected fainter host, with > 25.5 AB mag, 
may be located closer to the GRB position. 

To investigate the offset distribution in greater detail we 
supplement the values measured here with offsets for GRBs 
070724, 071227, and 0905 10 from g round-based observations 
dBergeret all 120091: iRauet al.ll2009l) . In the case of GRBs 
070724 and 071227 the optical afterglow s coincide with the 
disks of apparent edg e-on spiral galaxies dBerger et al.ll2009T: 
iD'Avanzo et al.ll2009l) . The offs ets of the three bursts are 4.8, 
14.8, and 5.5 kpc, respectively (iBerger et alj|2009t iRau et al.1 
120091) . For GRB 071227 we calculate the relative offset from 
our Magellan/IMACS observations and find a total (<70,grb + 
cre iga i) uncertainty of 65 mas, corresponding to 0.34 kpc at the 
redshift of the host 12 . 

There are 7 additional events with optical afterglow iden- 
tifications. Of these bursts, two (070707 and 070714b) co- 
incide with galaxies (Piranomonte et al. 2008; Graha met al.1 
2009), but their offsets have not been measured by the re- 
spective authors. Based on the claimed coincidence we con- 
servatively estimate an offset of < 0.5", corresponding 13 to 
< 4 kpc. Two additional bursts (070809 and 080503) do not 
have coincident host galaxies to deep limits, but the nearest 
galaxies are located abou t 6.5 and 20 kpc from th e afterglow 
positions, respectively 14 dPerley et a l. 2008, 2009). For the fi- 
nal three bursts (080905, 090305, and 090426) no deep host 
galaxy searches exist in the literature. 

In addition to the bursts with sub-arcsecond positions, 
several hosts have been identified within XRT error circles 
in follow-up observations (GRBs 060801, 061210, 061217 
070429b, 070729, and 080123; iBergeret all 120071: iBergeri 
2009), but in all of these cases the offsets are consistent 

12 This is signific antly more p recise than the large uncertainty of 0.4" 
quoted by D'Avanzo et al. 1 2009) based on absolute astrometry; for offset 
measurements differential astrometry provides a better approach. 

13 GRB 070714b is located at z = 0.923, while the redshift of GRB 070707 
is not known. Based on the faintness of the host, R ss 27.3 mag, we assume 
z = 1 to calculate the physical offset. 

14 GRB 070809 is located 19.6 kpc from a galaxy at z = 0.219, and about 
2.3" from a much fainter galaxy, which at z > 1 corresponds to 18.4 kpc. No 
host is detected at the position of GRB 080503 in deep HST observations, but 
a faint galaxy is located about 0.8" away, which at z > 1 corresponds to 6.5 
kpc. 



with z ero, or may be as large as ~ 30 kpc (e.g., Berg er et"ai1 
I2007T) . For example, the offsets for GRBs 060801, 061210, 
and 070429b are 19 ± 16 kpc, 11 ± 10 kpc, and 40 ±48 kpc. 
We use 30 kpc as a typical upper limit on the offset for these 
6 events. We note that no follow-up observations are avail- 
able in the literature for most short GRBs with X-ray positions 
from 2008-2009. Finally, about 1/4- 1/3 of all short GRBs 
discovered to date have only been detected in 7-rays, preclud- 
ing a unique host galaxy association and an offset measure- 
ment. 

The cumulative distribution of projected physical offsets 
for the GRBs with HST observations from this work, sup- 
plemented by the bursts with offsets or limits based on op- 
tical afterglow positions (070707, 070714b, 070724, 070809, 
071227, 080503, and 090510) is shown in Figure[l3] Also 
shown is the differential probability distribution, P{Sr)d{Sr), 
taking into account the non-Gaussian errors on the rad ial off- 
sets (see discussion in Appendix B of Bloom et al. 2002). We 
find that the median for this sample is about 5 kpc. 

As evident from the discussion above, this is not a complete 
offset distribution; roughly an equal number of short GRBs 
have only limits or undetermined offsets due to their detec- 
tion in just the X-rays or 7-rays 15 . Taking these events into 
account, our most robust inferences about the offset distribu- 
tion of short GRBs are as follows: 

• At least 25% of all short GRBs have projected physical 
offsets of < 10 kpc. 

• At least 5% of all short GRBs have projected physical 
offsets of > 20 kpc. 

• At least 50% of all short GRBs have projected physical 
offsets of < 30 kpc; this value includes the upper limits 
for the hosts identified within XRT error circles. 

These robust constraints are shown in Figure [T3l 

Using the observed distribution and the robust constraints 
outlined above, we now provide a comparison with predicted 
distributions for NS-NS binaries in Milky Way type galaxie s 
dBloom et al.ll999t iFrveret aljl999t iBelczvnski etal1l2006l) . 
appropriat e for the obs erved luminosities of short GRB host 
galaxies (Berger 2009|). We find good agree ment between 
the ob serve d distribution and those predicted by Bloo m et"aTI 
d 19991) and IBelczvnski et al.1 (2006). The offset distribution 
of iFrver et al.l d!999l) . with a median of about 7 kpc, predicts 
larger offsets and therefore provides a poorer fit to the ob- 
served distribution, which has a median of about 5 kpc. How- 
ever, all three predicted distributions accommodate the offset 
constraints. In particular, they predict about 60-75% of the 
offsets to be < 10 kpc, about 80-90% to be < 30 kpc, and 
about 10-25% of the offsets to be > 20 kpc. Thus, the pro- 
jected physical offsets of short GRBs are consistent with pop- 
ulation synthesis predictions for NS-NS binaries. However, 
the observations are also consistent with partial contribution 
from other progenitor systems with no expected progenitor 
kicks, such as WD-WD binaries. 

4.1. Host-Normalized Offsets 

To compare the offsets in a more uniform manner, we nor- 
malize the measured values by r e for each host galaxy. We 

15 We do not consider the bursts that lack host searches since there is no a 
priori reason that these events (mainly from 2008-2009) should have a differ- 
ent offset distribution compared to the existing sample from 2005-2007. 
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use the r e values measured from the one-dimensional radial 
surface brightness profiles from ellipse (see Figure fl4l 
and Table 01 and find values rangi ng from about 0.2 r e for 
GRB 060121 to 6.7 ± 2.7 r e for t helEvanset all d2009t) XRT 
position of GRB 051210. The iButled (120071) position for 
GRB 051210, however, leads to an offset of 4.65 ±4.60 r e , 
consistent with a negligible offset. For the subset of 6 bursts 
with optical afterglow positions and secure host associations, 
4 are located within 1 r e , while the remaining 2 bursts are 
located at about 2 r e (Figure [Hi. GRB 050509b, which has 
the largest physical offset, has a normalized offset of 2 - 3 r e , 
depending on which XRT position is used. Thus, with the ex- 
ception of the ambiguous case of GRB 061201, we find that 
all of the available offsets are consistent with < 2 r e . The 
large additional sample of physical offsets that we used above 
cannot be easily translated to host-normalized offsets at the 
present since none of the hosts have been observed with HST, 
thereby precluding a robust morphological analysis. This pro- 
vides an impetus for future HST observations. 

The differential probability distribution of host-normalized 
offsets for our HST sample, taking into account the non- 
Gaussian errors, is shown in Figure [15] We find that the me- 
dian value for all 8 bursts is ~ 1 r e . Moreover, < 20% of the 
probability distribution is at large offsets of > 2.5 r e . 

4.2. Comparison to Long GRB Host Galaxies 

We compare our observed sho rt GRB offse t s with those of 
long GRBs from the sample of iBloom et al.l (120021) in Fig- 
ures [T3~lfT5l The offset distribution of long GRBs has been 
used to argue for a mass ive star progenitor population, and 
against NS-NS binaries (IBloom et alJl2002t> . The offset dis- 
tribution for short GRBs is clearly shifted to larger physical 
scales. In particular, the median offset for the long GRBs is 
1.1 kpc, about a factor of 5 times smaller than the median 
value for short GRBs. Similarly, no long GRB offsets are 
larger than about 7 kpc, whereas at least some short GRBs 
appear to have offsets in excess of 15 kpc. 

However, the distinction between the two offset distribu- 
tions is significantly reduced when we take into account the 
systematically larger sizes of short GRB host galaxies (Fig- 
ure IT2"li. The median normalized offset for long GRBs is about 
0.8 r e , compared to about 1 r e for short GRBs. Similarly, 20% 
of the long GRB cumulative distribution has offsets of > 2.5 
r e , identical to the statistics for the short GRB offsets. Indeed, 
as can be see from Figure[l5] the cumulative host-normalized 
offset distributions for long GRBs and short GRBs with HST 
observations are nearly identical. 

In the context of NS-NS binary progenitors, the close simi- 
larity in the normalized offset distributions can be interpreted 
to mean that most systems likely remain bound to their hosts 
(rather than ejected into the intergalactic medium), and/or 
have a relatively short delay time. These conclusions are ten- 
tative due to the small number of events with host-normalized 
offsets, but they can be further tested with future HST obser- 
vations. 

5. LIGHT DISTRIBUTION ANALYSIS 

In addition to the offset analysis in the previous section, we 
study the local environments of short GRBs using a compari- 
son of their local brightness to the host light distribution. This 
approach is advantageous because it is independent of galaxy 
morphology, and does not suf fer from ambiguity in the def- 
inition of the host center (see iFruchter et al.1l2006l) . We note 
that for the overall regular morphology of short GRB hosts 



the definition of the h ost center is generally robust, unlike in 
the ca se of long GRBs (Frucht er et al.l20 06: Wainwrig ht et al.1 
120071) . On the other hand, this approach has the downside that 
it requires precise pixel-scale positional accuracy. In our sam- 
ple, this is the case for only 6 short bursts. 

The fraction of host light in pixels fainter than the afterglow 
pixel brightness for each host/filter combination is summa- 
rized in Table 2] The cumulative light distribution histogram 
is shown in Figure [16] The shaded histogram represents the 
range defined by the dual filters for 5 of the 6 bursts. We find 
that the upper bound of the distribution is defined by the blue 
filters, indicating that short GRBs trace the rest-frame opti- 
cal light of their hosts better than the rest-frame ultraviolet. 
This indicates that short GRB progenitors are likely to be as- 
sociated with a relatively old stellar population, rather than a 
young and UV bright population. 

The overall distribution has a median value of ss 0.1 -0.4 
(red); namely, only in about one-quarter of the cases, 50% 
of the host light is in pixels fainter than at the GRB location. 
Thus, the overall distribution of short GRB locations under- 
represents the host galaxies' light distribution. This is also 
true in comparison to the distribution f or core-collapse SNe , 
which appear to track their host light (Fruc hter et al.l 2006). 
and even Type la SNe, which have a median of about 0.4 
dKellv et alj|2008l) . Thus, the progenitors of short GRBs ap- 
pear to be more diffusely distributed than Type la SN progen- 
itors. 

5.1. Comparison to Long GRB Host Galaxies 

An extensive analysis of the brightness distrib ution at the 
locatio n of long GRBs has been carried out by IFruchter et al.1 
(2006). These authors find that long GRBs are more con- 
centrated on the brightest regions of their hosts than expected 
from the light distribution of each host. In particular, they 
conclude that the probability distribution of GRB positions is 
roughly proportional to the surface brightness squared. As 
can be seen from Figure [16] short GRBs have a significantly 
more diffuse distribution relative to the host light than long 
GRBs. In particular, for the latter, the median light fraction is 
about 0.85 compared to about 0.25 ±0.15 for the short GRBs. 

6. DISCUSSION AND IMPLICATIONS 

Our extensive analysis of short GRB host galaxy morpholo- 
gies and the burst local environments has important implica- 
tions for the progenitor population. We address in particular 
the popular NS-NS merger model, as well as delayed magne- 
tar formati on via WD- WD me rgers or WD accretion-induced 
collapse dMetzger et al.ll2008l) . 

6.1. Morphology 

From the morphological analysis we find continued evi- 
dence that the bulk of short GRB host galaxies (~ 3/4) are 
late-type galaxies, in agreement w ith results from spectro- 
scopic observations (Berger 2009). Moreover, as demon- 
strated by the systematic differences in luminosity, star for- 
mation rates, and metalli cities betwee n the star forming hosts 
of long and short GRBs (Berger 2009), we find here that short 
GRB hosts are systematically larger than long GRB hosts. 
These results indicate that the progenitors of the two GRB 
classes select different environments. The higher luminosi- 
ties, larger sizes, and lower specific star formation rates of 
short GRB hosts suggest that their rate of occurrence is tied 
to galactic mass rather than to star formation activity. This 
result is in broad agreement with old progenitor populations 
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such as NS-NS, NS-BH, or WD-WD binaries, but it indicates 
that the bulk of short GRB progenitors are not young magne- 
tars. This conclusion is also supported by the dearth of merger 
signatures, which point to delays of > 10 9 yr relative to any 
merger-triggered star formation episodes. 

6.2. Offsets 

The differential offsets measured here from the HST ob- 
servations provide the most precise values to date for short 
GRBs, with a total uncertainty of only ~ 10-60 mas, corre- 
sponding to ~ 30-500 pc. We find that none of the offsets 
are smaller than ~ 1 kpc, while this is the median offset for 
long GRBs. On the other hand, a substantial fraction of the 
measured offsets are only a few kpc. The median offset for 
the HST observations supplemented by ground-based data is 
about 5 kpc (Figure [T3l . about 5 times larger than for long 
GRBs. 

As discussed in detail in $4] the observed offset distribution 
is incomplete. About 1/4-1/3 of all short GRBs have only 
7-ray positions (~ 1 — 3'), and a similar fraction have only 
XRT positions, which generally lead to a range of offsets of 
^0-30 kpc. Taking these limitations into account we find 
that the most robust constraints on the offset distribution are 
that > 25% of all short GRBs have offsets of < 10 kpc, and 
that > 5% have offsets of > 20 kpc. Both the observed off- 
set distribution and these constraints are in good agreement 
with predictions for the offset distribution of NS-NS binaries 
in Milky Way type galaxies dBloom et al.l Il999t iFrver et all 
119991; [Belczvns ki et al.l 12006b . However, at the present they 
cannot rule out at least a partial contribution from other pro- 
genitor systems such as delayed magnetar formation and even 
young magnetar flare. The apparent existence of large offsets 
in the sample suggests that these latter models are not likely 
to account for all short GRBs. 

In contrast to the larger physical offsets of short GRBs, 
we find that the distribution of offsets normalized to the host 
galaxy effective radii exhibits much better agreement between 
long and short GRBs (Figure [15). The medians of the two 
distributions are similar (1 versus 0.8 r e for short and long 
GRBs, respectively), and both populations have k, 20% prob- 
ability for offsets of > 2.5 r e . Naturally, due to the lack of 
HST observations for short GRBs from 2007-2009, the sam- 
ple for which this analysis is possible is smaller than the sam- 
ple with physical offsets. Thus, further HST observations of 
existing and future short GRB hosts are essential in order 
to determine whether the broad similarity in host-normalized 
offsets is robust. We stress that in the context of comparing 
short GRBs with various progenitor populations and with long 
GRBs, host-normalized offsets are the more relevant quantity. 
For example, the host-normalized distributions for massive 
stars in small and large galaxies will be similar even though 
their physical offsets will differ. We stress that population 
synthesis modelers should include an appropriate distribution 
of host galaxy sizes, and thereby provide predictions for host- 
normalized offsets. 

In the context of implications for the progenitor population, 
a rece nt study of short GRB physical offsets by iTroia et al.l 
(2008) led these authors to claim that short GRBs with ex- 
tended X-ray emission have systematically smaller offsets, 
possibly due to a systematic difference in the progenitors. Our 
HST sample includes three short GRBs with strong extended 
emission (050709, 050724, and 061006), and one GRB 
(060 121) with possib le extended emission (4.5cr significance; 
Donaghv et al. 2006). The physical offsets of these bursts are 



about 3.7, 2.7, 1.3, and 1 kpc, respectively, leading to a mean 
offset of about 2.2 kpc. The physical offsets of the bursts with- 
out extended emission, but with precise afterglow positions 
(051221, 060313, and 061201) are 2.0, 2.3, and 14.2 or 32.5 
kpc, respectively. The two events with no extended emission 
and with XRT positions (050509b and 051210) have offsets 
of about 54 ± 12 and 28 ± 23 kpc, respectively. If we include 
the ground-based sample with optical afterglow positions (see 
$4), we find that the bursts with apparen t extended emission 
(070714b, 071227, 080513, and 090510;|Barbier et al. l )2007) ; 
ISakamoto et alj|2007l;lukwatta et alj2009l;lPerlev et alj2009i) 
have offsets of < 4, 14.8, ~ 20, and ^5.5 kpc, while the 
bursts without extended emission (070724 and 070809) have 
offsets of 4.8 and ~ 6.5 kpc. Thus, based on the sample of 
events with sub-arcsecond positions we find that 6/8 bursts 
with extended emission have offsets of < 5 kpc and 2/8 have 
likely offsets of ~ 15-20 kpc. In the sample without extended 
emission we find that 4/5 have offsets of < 6 kpc and 1/5 has 
a likely offset of ~ 14-32 kpc. Thus, we conclude that there 
is no significant difference in the two offset distributions. 

The inclusion of events with only XRT positions does not 
change this conclusion. In particular, of the subset with no 
extended emission only GRB 050509b is likely to have a sig- 
nificant offset, while GRBs 051210, 060801, and 070429b 
have offsets (28 ± 23, 19 ± 16, and 40 ± 48 kpc, respectively) 
that are consistent with zero. Similarly, GRB 061210 with ex- 
tended emission has an offset of 1 1 ± 10 kpc. An examination 
of the sample of ITroia et alj d2008l) reveals that their claim 
that short GRBs without extended emission have systemati- 
cally larger offsets rests on four events in particular: GRBs 
050509b, 060502b, 061217, and 061201. As noted above, 
GRBs 050509b and 061201 indeed appear to have substan- 
tial offsets 16 , but so do GRBs 071227 and 080503 with ex- 
tended emission and offsets of about 15-20 kpc. Next, the 
large offset for GRB 060502b relies on its claimed associa- 
tion with an elliptical ga laxy 70 ± 16 kpc from the XRT posi- 
tion dBloom et al.ll200 7). However, the XRT erro r circle con- 
tains additional galaxies with negligible offsets (Berg er et alj 
12007b . Finally, we note that the offset for GRB 061217 is un- 
reliable due to a sub stantial discr epan cy of about 33" in the 
XRT positions from lButleri (120071) and lEvans etail (l2009t) . A 
continued investigation of the difference between short GRBs 
with and without extended emission will greatly benefit from 
the use of host-normalized offsets. 

6.3. Light Distribution 

In addition to projected offsets relative to the host center, 
we find that the locations of the short GRBs with HST imag- 
ing and sub-arcsecond positions are more diffusely distributed 
relative to their host light than long GRBs. In particular, we 
find that short GRB positions under-represent their host light, 
even in comparison to core-collapse and Type la SNe. This 
result is likely an upper limit on the brightness of short GRB 
locations since only the subset of events with optical after- 
glow positions can be studied with this approach. Thus, short 
GRBs arise from a population of events with a more diffuse 
distribution than massive stars and Type la SN progenitors. 
This result also indicates that the bulk of the progenitors of 
long and short GRBs cannot both be magnetars. 

16 We note that lTroia et aH <200l) assume that the host of GRB 061201 is 
the galaxy at z = 0.111 at an offset of about 32 kpc. However, as we have 
shown here based on the HST observations (Appendix [5), there is a fainter 
potential host at a smaller offset. 
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There are currently 10 known short GRBs with optical af- 
terglows for which HST observations will enable a similar 
analysis. This is twice the number of the current sample, and 
we can therefore make significant progress in understanding 
the relation of short GRB environments to the overall distri- 
bution of light in their host galaxies with future observations. 

7. CONCLUSIONS 

We presented the first comprehensive analysis of short GRB 
HST observations, and used these data to extract the morpho- 
logical properties of the host galaxies, the projected physical 
and host-normalized GRB offsets, and the brightness at the 
location of the bursts relative to the overall light distribution 
of their hosts. The main conclusions of our analysis are as 
follows: 

• The majority of short GRB hosts are consistent with or 
have exponential surface brightness profiles, typical of 
late-type galaxies. This conclusion is in good agree- 
ment with results from spectroscopic observations that 
reveal star formatio n activity in ~ 3/4 of short GRB 
hosts (lBergeril2009l) . 

• The host galaxies of short GRBs are on average larger 
by about a factor of 2 than the hosts of long GRBs. 

• The observed short GRB offset distribution extends 
from ~ 1 to 50 kpc, with a median of about 5 kpc. In- 
cluding the short GRBs with only 7-ray or X-ray posi- 
tions, we find that > 25% of all events have offsets of 

< 10 kpc, and > 5% have offsets of > 20 kpc. A addi- 
tional, though softer, limit is that > 50% have offsets of 

< 30 kpc. 

• The observed physical offset distribution and the robust 
constraints compare favorably with the predicted distri- 
bution for NS-NS binaries. However, they do not rule 
out at least a partial contribution from other progenitors 
systems such as WD-WD binaries. 



• We find no convincing evidence that short GRBs with 
extended emission have smaller physical offsets than 
those without extended emission. In both sub-samples 
we find examples of both small offsets (~ few kpc) and 
possibly large offsets (tens of kpc). 

• The distribution of host-normalized offsets for the sub- 
set of short GRBs with HST observations is nearly iden- 
tical to that of long GRBs. This is due to the system- 
atically larger size of short GRB hosts, and indicates 
that a comparison with long GRBs and progenitor mod- 
els will benefit from the use of host-normalized (rather 
than physical) offsets. 

• The locations of short GRBs with sub-arcsecond posi- 
tions and HST imaging under-represent the overall light 
distribution of their hosts, but less so in the red. This re- 
sult differs substantially from long GRBs, core-collapse 
SNe, and even Type la SNe. 

The results derived in this paper are based mainly on a small 
sample of short GRBs (9 events) from 2005-2006. Seven of 
these objects have precise positions based on optical after- 
glow detections. Ten additional events with precise afterglow 
positions, and a similar number with XRT positions (some of 
which with identified hosts), are now available for a similar 
study. It is essential to observe this existing sample, as well 
as new events from Swift and Fermi, with the refurbished HST 
using the ACS and WFPC3 instruments. In conjunction with 
constraints on the progenitor population from the redshift dis- 
tribution (Berger et al. 2007) and spectroscopic studies of the 
host galaxies (lBerger]2009t) . the continued use of high angular 
resolution imaging will provide crucial insight into the nature 
of the progenitors and the potential for multiple populations. 



APPENDIX 
THE HOST GALAXY OF GRB 0603 13 

We present the first host galaxy association for GRB 060313, using HST/ ACS observ ations in the F47 5W and F775W filters. 
The offset between the GRB position (determined from Gemini-South observations; Berger et al. 2007) and the galaxy center 
is about 0.32" (Table |3). The galaxy brightness is m(F475W) = 26.4 AB mag and m(F775W) = 25.6 AB mag (Tab le [p. The 
probability of chance coincidence at this offset and galaxy brightness is only about 3 x 10~ 3 ((Beckwith et al. 2006). We thus 
conclude that this galaxy is the likely host of GRB 060313. 

POSSIBLE HOST GALAXIES OF GRB 061201 

The HST observations of GRB 061201 and its environment are shown in Figure [10] We explore two possibilities for t he host 
galaxy. First, the b urst is located 16.2" (32.5 kpc) from a relatively bright galaxy at z = 0. 1 1 1 (marked "A" in Figure[l0j Berger 
2006; Stratta et alj|2007l) . for which we measure m(F606W) = 18.17 and m(F814W) = 17.82 AB mag. Second, we identify from 
the HST/ ACS observations a second, fainter galaxy (marked "B" in Figure [TUb located 1.8" from the GRB position, and with 
ra(F606W) = 25.34 and ra(F814W) = 25.03 AB mag. The redshift of this galaxy is not known, but as suming z > 1 the inf erred 
projected offset is 14.2 kpc. The probability of chance coincidence for both galaxies is about 20% dBeckwith et al J 120061) . We 
therefore do not claim a unique host galaxy association for this burst, and stress that both galaxies should be considered as 
potential hosts. Deeper HST observations may also uncover an underlying host. 

POSSIBLE HOST GALAXIES OF GRB 060502B 

The HST observations of GRB 060502b and its environment are shown in Figure [8] Previously. iBloom et alj (120071) claimed 
that the host is likely an early type galaxy at z = 0.287 located about 70 kpc away from the burst XRT position. These authors 
also no te the presence of fainter objects within the XRT error circle. A galaxy with R ks 25.2 mag was also found by Berge r et"aTI 
d2007l) . In the combined //S77ACS/F814W we find 6 faint galaxies within the XRT error circles of GRB 060502b (Figure [8]). 
These galaxies have the following AB magnitudes: 27.5 (B), 25.9 (C), 27.2 (D), 26.1 (E), 24.8 (F), 25.7 (G). The probability of 
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chance coincidence for these galaxies within the XRT error circles is of the order of unity. 
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TABLE 1 

HST Observations of Short GRB Host Galaxies 



GRB 


RA 
(J2000) 


Dec 
(J2000) 


Uncert. 

(") 


OA? 




Instrument 


Filter 


Date 
(UT) 


Exp. Time 

(s) 


AB mag fl 


4 h 

(mag) 


050509b 


12 h 36 m 14.06 s 
12 h 36 m 13.76 s 


+28°59'07.2" c 
+28°59'03.2" 


3.4 
3.3 


N 


0.226 


ACS 


F814W 


2005 May 14 


6870 


16.32 


0.037 


050709 




JO JO Jy.J 


2 


Y 


1606 


ACS 
WFPC2 


F814W 
F450W 


9006 Till \fy 

2007 Jul 29 


6981 
3200 


21 09 
21.43 


02 
0.045 


050724 


16 h 24 m 44.38 s 


-27°32'27.5" 


0.1 


Y 


0.257 


WFPC2 
WFPC2 


F450W 
F814W 


2008 Apr 07 
2008 May 18 


3200 
3200 


19.98 
18.74 


2.645 
1.189 


051210 


22 h 00 m 41.26 s 
22 h 00 m 41.33 s 


-57°36'46.5" 
-57°36'49.4" 


2.9 
1.7 


N 


> 1.4 


WFPC2 


F675W 


2007 Apr 03 


2800 


21.14 


0.052 


051221a 


21 h 54 m 48.63 s 


+16°53'27.4" 


0.2 


Y 


0.5465 


WFPC2 
WFPC2 


F555W 
F814W 


2007 Aug 13 
2007 Aug 22 


3200 
1600 


21.86 
21.42 


0.227 
0.133 


060121 


09 h 09 m 51.99 s 


+45°39'45.6" 


0.1 


Y 




ACS 


F606W 


2006 Feb 27 


4400 


26.22 


0.047 


060313 


04 h 26 m 28.42 s 


-10°50'39.9" 


0.2 


Y 




ACS 
ACS 


F475W 
F775W 


2006 Oct 13 
2006 Oct 14 


2088 
2120 


26.38 
25.61 


0.300 
0.135 


060502b 


18 h 35 m 45.53 s 
18 h 35 m 45.28 s 


+52°37'52.9" 
+52°37'54.7" 


3.7 
5.8 


N 




ACS 


F814W 


2006 May 15-Jul 16 


25224 


17.88/ 24.8-27.5 d 


0.085 


061006 


07 h 24 m 07.78 s 


-79°11'55.5" 


0.2 


Y 


0.4377 


ACS 
WFPC2 


F814W 
F555W 


2006 Oct 14 
2008 May 22 


6054 
3200 


21.67 
23.90 


0.616 
1.052 


061201 


22 h 08 m 32.09 s 


-74°34'47.1" 


0.2 


Y 


0.111/ ••• e 


ACS 
ACS 


F606W 
F814W 


2006 Dec 1 1 
2006 Dec 11 


2178 
2244 


18.17 /25.34^ 
17.82/25.03'' 


0.251 
0.147 



NOTE. — Summary of short GRB positions and redshifts, HST observations, and host galaxy magnitudes (calculated using IRAF/ellipse). 
a These values have been corrected for Galactic extinction. 

b Galactic extinction. 

c In all cases with Swift/XRT positions the top and bottom set of coordinates are from the catalogs of Butler 1 2007) and Evans et al. ( 2009), respectively. 

These magnitudes correspond to galaxy "A" and galaxies '*B*'- "G" in Figure [8l 
e We consider two possible host galaxies for this burst (see Appendix^}. 
* The first value is for galaxy "A" and the second is for galaxy "B" in Figure fTO] 
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TABLE 2 

Morphological Properties of Short GRB Host Galaxies 



galfit IRAF/ellipse 



GRB 


Instrument 


Filter 


n" 


r e 




h 

Me 


AB Mag* 


n 


I'e 














(") 


(kpc) 


(AB mag arcsec 






(") 


(kpc) 


(AB mag arcsec 


U J UD Uy D 




CO 1 A\U 

rol'tw 


S ft 






11 S 


i f* i 


J.O 


J.Ot 


zu.yo 


11 A 


050709 


ACS 


F814W 


1.1 


0.76 


2.08 


23.0 


21.2 


0.6 


0.64 


1.75 


22.4 




W rrLZ 


r4juw 


1 1 


u. / 1 


1 Q -1 


11 ^ 


1 Q 

zi.y 


ft Q 

u.y 


n 7 1 

U. / 1 


1 Q -1 


11 A 


050724 


WFPC2 


F450W 


4 


1.35 


5.34 


23.2 


19.4 


1.3 


0.36 


1.42 


20.8 




WFPC2 


F814W 


3.0 


0.82 


3.24 


20.5 


18.0 


2.9 


1.01 


4.00 


20.8 


05 1 2 1 


WFPC2 


r\j i j w 


j 


70 


J.Oj 


23 8 


91 7 


1 


U.U J 


5 07 


24 2 








4 


2.38 


19.14 


26.3 


22.8 










051221a 


WFPC2 


F555W 


0.9 


0.36 


2.29 


23.3 


23.1 


0.8 


0.34 


2.17 


23.1 




WFPC2 


F814W 


0.9 


0.41 


2.61 


22.7 


22.1 


0.9 


0.39 


2.49 


22.7 


060121 


ACS 


F606W 


1 


0.36 


2.89 


25.9 


27.1 


1.4 


0.67 


5.39 


27.2 








4 


1.22 


9.81 


27.4 


26.6 










060313 


ACS 


F475W 


1 


0.14 


1.13 


23.7 


27.3 


0.6 


0.17 


1.37 


24.9 








4 


0.32 


2.57 


26.2 


26.7 












ACS 


F775W 


1 


0.07 


0.56 


21.4 


26.1 


1.3 


0.23 


1.85 


25.0 








4 


0.10 


0.80 


23.6 


26.3 










061006 


ACS 


F814W 


0.7 


0.57 


3.22 


22.3 


22.7 


0.7 


0.65 


3.67 


22.9 




WFPC2 


F555W 


1 


0.63 


3.55 


23.3 


23.4 


0.8 


0.55 


3.10 


23.6 



NOTE. — Results of morphological analysis performed with galfit and IRAF/ellipse (g3|. 
a For the Sersic index n in galfit, exact values of 1 and 4 indicate a fit with n as a fixed parameter. 
b These values have been corrected for Galactic extinction. 

c For the hosts with unknown redshift (GRBs 051210, 060121, and 060313) we assume z = I. 
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TABLE 3 

Short GRB Angular, Physical, and Host-Normalized Offsets 



GRB 


Instrument 


Filter 


z 


Reference 


No. 


^GB^HST 


°"e,GRB 


°0,gal 


<5RA 


<5Dec 


Offset 


Offset 


Offset 














(mas) 


(mas) 


(mas) 


(") 


(") 


(") 


(kpc) 


(r e °) 


050509h 


ACS 


F814W 


0.226 


SDSS 


10 


30 


3400 


1 .0 


+15.61 


+8.40 


17.73 ± 3.4 


63 7 _j_ 12.2 


3 04 + 58 

. \J*T _1_ U.JO 
















3300 


1.0 


+11.68 


+4.40 


12.48 ±3.3 


44.8 ± 11.9 


2.14±0.57 




Arc 


roI4 W 


n i aha 

U. I DUO 




3^ 

j j 


Q 

o 


i n 
1 ,U 


1 A b 
1.4 




a Tin 


1 111 | n nin 
I.jjI ± U.U1U 


J.D4 it U.Uz / 


o no _l n m 




WFPC2 


F450W 




HST/ACS 


12 


14 


1.0 


4.7 


+1.306 


-0.360 


1.355 ±0.020 


3.71 ±0.055 


1.91 ±0.03 




ACS 


F814W 




self 






1.0 


1.4 


+1.329 


-0.310 


1.365 ±0.002 


3.74 ±0.005 


2.13±0.01 


050724 


WFPC2 


F450W 


0.257 


Magellan/PANIC 


60 


15 


5.0 


4.7 c 


-0.226 


-0.640 


0.679 ±0.025 


2.69 ±0.099 


1.89 ±0.07 




WFPC2 


F814W 




Magellan/PANIC 


85 


14 


5.0 


1.4 


-0.213 


-0.630 


0.665 ±0.020 


2.63 ±0.079 


0.66 ±0.02 




PANIC 


K s 




self 






5.0 


1.0 


-0.253 


-0.650 


0.697 ±0.006 


2.76 ±0.024 




051210 


WFPC2 


F675W 


> 1.4 


2MASS 


12 


29 


2900 


8 


+2.89 


+0.50 


2.93 ±2.9 


24.9 ± 24.6 


4.65 ±4.60 
















1700 


8 


+3.45 


-2.40 


4.20 ± 1.7 


35.7 ± 14.4 


6.67 ±2.70 


051221 


WFPC2 


F555W 


0.5465 


Gemini-N/GMOS 


45 


23 


2.5 


3.1 


-0.287 


+0.090 


0.301 ±0.029 


1.92±0.18 


0.88 ±0.08 




WFPC2 


F814W 




Gemini-N/GMOS 


45 


24 


2.5 


3.1 


-0.330 


+0.090 


0.342 ±0.030 


2.18±0.19 


0.88 ±0.08 


060121 


ACS 


F606W 




Gemini-N/GMOS 


25 


18 


16 


12 


-0.115 


+0.030 


0.1 19 ±0.046 


0.96±0.37 < ' 


0.18 ±0.07 


060313 


ACS 


F475W 




Gemini-S/GMOS 


30 


30 


19 


19.4 


+0.354 


+0.040 


0.356 ±0.068 


2.86±0.55'' 


2.09 ±0.40 




ACS 


F775W 




Gemini-S/GMOS 


15 


30 


19 


13.2 


+0.280 


+0.050 


0.284 ±0.062 


2.28 ±0.50'' 


1.23 ±0.23 


060502b 


ACS 


F814W 




USNO-B 


47 


120 


3700 




























5800 














061006 


ACS 


F814W 


0.4377 


VLT/FORS1 


75 


17 


21 


4.9 


-0.155 


-0.170 


0.230 ±0.043 


1.30 ±0.24 


0.35 ±0.07 




WFPC2 


F555W 




VLT/FORS1 


45 


20 


21 


11 


-0.171 


-0.190 


0.256 ±0.052 


1.44 ±0.29 


0.46 ±0.10 


061201 


ACS 


F814W 




VLT/FORS2 


24 


13 


41 
















ACS 


F606W 




VLT7FORS2 


24 


13 


41 















NOTE. — Projected angular, physical, and host-normalized offsets for the short GRBs with HST observations. 
" Values for r e are from ellipse (Table[2l 
h Systematic uncertainty in host center is 50 mas. 
c Systematic uncertainty in host center is 20 mas. 
d Assuming z= 1. 

e We do not claim a unique host galaxy identification for this burst. 



HST Observations of Short GRB Host Galaxies 



TABLE 4 
Short GRB Fractional Flux 



GRB 


Instrument 


Filter 


Fractional Flux 




WFPC2 


rnju vv 


o 




ACS 


F814W 


0.09 


0S0794. 


WFPP9 
vv m,z 


r^ju vv 


n rn 

U.UJ 




WFPC2 


F814W 


0.33 


051221 


WFPC2 


F555W 


0.54 




WFPC2 


F814W 


0.65 


060121 


ACS 


F606W 


0.41 


060313 


ACS 


F475W 


0.04 




ACS 


F775W 





061006 


WFPC2 


F555W 


0.56 




ACS 


F814W 


0.63 



NOTE. — Fraction of host galaxy light in pixels fainter 
than the GRB position. 
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GRB050509b 




Fig 1 .— Le ft: HSTI ACS /F814W image of the location of GRB 050509b. The circles mark the X-ray positions of the afterglow from the analysis of lButleil 
12007) (red) and Evans et al. (2009) (blue). Center: Sersic model fit from galf it. Right: Residual image. 
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N 

J 



2.74 kpc 



FIG. 2. — Top-left: //S7VACS/F450W image of the host galaxy of GRB 050709 with a 10cr error circle representing the afterglow position. Top-center: Sersic 
model fit from galf it. Top-right: Residual image. Bottom: Same, but for the //S77WFPC2/F814W observations. 
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FIG. 3 . — Top-left: ffi>r/WFPC2/F450W image of the host galaxy of GRB 050724 with a 5cr error circle representing the afterglow position. Top-center: Sersic 
model fit from galf it. Top-right: Residual image. Bottom: Same, but for the //S77WFPC2/F814W observations. 



HST Observations of Short GRB Host Galaxies 1 7 




FIG. 4. — Top-left: //S77WFPC2/F675W image of the location of GRB 051210. The circles mark the X-ray positions of the afterglow from the analysis of 
IButleJ (2007) (red) and lEvans et al. (2009) (blue). Top-center: Sersic model fit from gal fit with a fixed value of n = 1. Top-right: Residual image. Bottom: 
Same, but for Sersic model fit from galf it with a fixed value of n = 4. 
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GRB051221a 








WFPC2/F555W 








MY. 


• 




N 

,J 


i" 

6.37 kpc 








GRB051221a 






mm Qj 


WFPC2/F814W 
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_ 




N 

,J 


i" 

6.37 kpc 







FIG. 5. — Top-left: ffi>r/WFPC2/F555W image of the host galaxy of GRB 051221 with a 5cr error circle representing the afterglow position. Top-center: Sersic 
model fit from galf it. Top-right: Residual image. Bottom: Same, but for the //S77WFPC2/F814W observations. 
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GRB060121 




Mi 

GRB060121 
VDS/F606W, 




■I- 1 5 1 ■ rf-" jl 




FIG. 6. — Top-left: HST/ ACS/F606W image of the host galaxy of GRB 060121 with a 3cr error circle representing the afterglow position. The image has been 
smoothed with a 2 X 2 pixel Gaussian. Top-center: Sersic model fit from galf it with a fixed value of n = 1. Top-right: Residual image. Bottom: Same, but for 
Sersic model fit from galf it with a fixed value of n = 4. 
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GRB060313 
ACS/F475W 




GRB060313 
ACS/F775W 



FIG. 7. — Top-left: //S77ACS/F475W image of the host galaxy of GRB 060313 with a 3cr error circle representing the afterglow position. The image has been 
smoothed with a 3 X 3 pixel Gaussian. Top-center: Sersic model fit from galf it. with a fixed value of n = 1, Top-right: Residual image. Second panel: Same, 
but for Sersic model fit from galf it with a fixed value of n - 4. Third panel: Same, but for the //S77WFPC2/F775W observations with a fixed value of n = 1. 
Bottom panel: Same, but for the //S77WFPC2/F775W observations with a fixed value of n = 4. 
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FIG 8. — g,S T/ACS/F814W imag es of the host galaxy of GRB 060502b. The circles mark the X-ray positions of the afterglow from the analysis of lButleri 
120071) (red) and lEvans et"aT] <2009t ) (blue). The bright galaxy marked "A" is located at z = 0.287 IBloom et af]|2007D . Several fainter galaxies ("B"-"G") are 
located within the XRT error circles (see Appendix O. 
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FIG. 9. — Top-left: //ST/ACS/F555W image of the host galaxy of GRB 061006 with a 3<r error circle representing the afterglow position. Top-center: Sersic 
model fit from galf it. Top-right: Residual image. Bottom: Same, but for the ffi>77ACS/F814W observations. 
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GRB061201 
ACS/F606W 




FIG . 10.— Left: HST/ A CS/F606W image of the location of GRB 06 1 20 1 . The bright galaxy at the upper right-hand corner ("A") is located at z = 0. 1 1 1 I Berger 
2006; Stratta et al. 2007) with an offset of about 32.5 kpc. A second, fainter galaxy ("B") is located about 1.8" away from the optical afterglow position (see 
Appendix|B). Right: Same, but for the //S77ACS/F814W observations. 
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Semi-major Axis (arcsec) 

FIG. 1 1 . — One-dimensional radial surface brightness profiles for short GRB host galaxies derived from IRAF/ellipse. The gray lines are Sersic model fits 
(Equation[T} to the surface brightness profiles. The results of the fits are listed in Tableff] 
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Sersic n Fraction 

FIG. 12. — Effecti ve ra dii for the short GRB hosts with HST observations plotted as a function of their Sersic n values. We use the results of the IRAF/ellipse 
analysis; see Figure [TT| (open squares designate hosts for which q alf it models with n = 1 and n = 4 provide an equally good fit). Also shown are the data 
for long GRB hosts based on HST observations from the sample of Wainwright et al. 1 2007). The hosts of GRBs 050509b and 050724 have n values typical of 
elliptical galaxies, but the remaining hosts have a similar distribution to that of long GRBs (i.e., a median of n ~ 1, or an exponential disk profile). On the other 
hand, the hosts of short GRBs are larger by about a factor of 2 than the hosts of short GRBs, in agreement with their higher luminosities. 
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Offset (kpc) 

FIG. 13. — Projected physical offsets for short GRBs (black) and long GRBs (gray; Bloom et al. 2002). The top panel shows a cumulative distribution, while 
the bottom panel shows the differential distribution taking into account the non-Gaussian errors on the offsets. The arrows in the bottom panel mark the median 
value for each distribution. The median value for short GRBs, ?s 5 kpc, is about a factor of 5 times larger than for long GRBs. The arrows in the top panel 
exhibit the most robust constraints on the offset distribution (i|4), taking into account the fraction of short GRBs with only 7-ray positions, as well as short GRBs 
for which hosts have been identified within XRT error circles (thereby providing a typical range of ~ — 30 kpc). Also shown in the top panel are predicted 
offset distributions for NS-NS binary mergers in Milky Way type galaxies based on population synthesis models. We find good agreement between the observed 
distribution and models, as well as between the robust constraints and models. 
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Short GRBs 
Long GRBs 




FIG. 14. — Host-normalized offsets in right ascension and declination for the short GRBs in our HST sample (black; open symbols mark the GRBs with X-ray 
positions, 050509b and 051210). Also shown are the offsets for long GRBs from the sample of Bloom et al. 12002). The circle marks an offset of 1 r e . About 
half of all long GRBs have offsets of < 1 r e , and we find a similar result for short GRBs. 
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T 




log (5r / g 

FIG. 15. — Probability distributions of the host-normalized offsets of short GRBs (black) and long GRBs (gray; B loom et al. 2002). For each burst we include 
the host-normalized offset taking into account the non-Gaussian errors. The arrows mark the median value of each distribution, and the inset shows the cumulative 
distribution. 
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FIG. 16. — Cumulative distribution of fractional flux at the location of short GRBs relative to their host light. For each burst we measure the fraction of host 
light in pixels fainter than the GRB pixel location. The shaded area is defined by the results for th e two available filter s for each s hort GRB. Also shown are data 
for long GRBs (dark gray line) and for core-collapse and Type la SNe (light gray lines) from Fruchter et al. 1 2006) and Kellv et al. 1 2008). The dashed line marks 
the expected distribution for objects which track their host li ght distribu tion. Short GRBs appear to under-represent their host light, while long GRBs tend to be 
concentrated in the brightest regions of their hosts (Fruchter et al. 2006). 



